Localized plasmons in point contacts 
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^ \ Abstract 

^ I Using a hydrodynamic model of the electron fluid in a point contact geometry we show that 
{Sj | localized plasmons are likely to exist near the constriction. We attempt to relate these plasmons 
with the recent experimental observation of deviations of the quantum point contact conduc- 
tance from ideal integer quantization. As a function of temperature this deviation exhibits an 
activated behavior, exp(— T a /T), with a density dependent activation temperature T a of the 
order of 2 K. We suggest that T a can be identified with the energy needed to excite localized 
1) | plasmons, and we discuss the conductance deviations in terms of a simple theoretical model 
involving quasiparticle lifetime broadening due to coupling to the localized plasmons. 
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Introduction 

The quantized conduction through a narrow point contact is one of the key effects in mesoscopic 
physics, the quantum point contact remains an important testing ground for the description 
O ' of mesoscopic phenomena. Recently, significant deviations from the Landauer-Biittiker theory 
i^i - have been observed in quantum point contacts in the temperature dependence of the conduc- 
. tance quantization [[I], 0] and as a so called "0.7" structure or quasi plateau, appearing around 
j> ! 0.7 times the conductance quantum 2e 2 /h ||. Invoking a Luttinger liquid approach 0] the de- 
CN \ viations have been discussed in terms of interaction effects []5], [| [/J. However, firm conclusions 
^ 1 have been difficult to obtain partly due to the narrow temperature range (0.1 K - 4 K) in which 
t^- ■ the effect can be studied in conventional split gate quantum point contacts, where relatively 
■ close lying one-dimensional subbands are formed. 

! An important progress was provided with the appearance of strongly confined GaAs quan- 
tum point contacts using a combination of shallow etching and a top gate ||. In these new 
samples the conduction quantization can be followed up to around 30 K. In a subsequent work 
H these samples were used to study the temperature dependence deviations from perfect con- 
ductance quantization. At low temperature (~ 0.05 K) almost ideal quantized conductance 



O . is observed for the first conduction plateau, but deviations develop as the temperature is in- 
*P. \ creased. The enlarged temperature range allowed for the observation of activated temperature 
| dependence of these deviations: 5G(T) oc exp(— T a /T). Furthermore, by changing the top gate 
^ ■ it was found that T a increases with increasing density. An explanation could not be found 
■ using the standard single particle picture, and in the brief theory section of Ref. || we there- 
fore suggested to include collective effects through plasmons. In short, we identified T a as the 
energy needed to excite localized plasmons, and we discussed the conductance in terms of a 
simple theoretical model involving the additional effect of electrons scattering off the localized 
plasmons. In the present theoretical work we elaborate on that idea. In a point contact the 
charge is of course depleted. In order to study the collective excitations of such a system, we 
can approach it from two limits: 1) starting from a homogeneous electron liquid which is 0, 
or 2) starting with two spatially separated liquids. Below we follow the first approach, and we 
argue from a hydrodynamic model that localized plasmons may exist in realistic situations. 



Plasmons of a homogeneous electron liquid in a cylinder 

Following Fetter [K| we use a hydrodynamic model of a weakly damped, compressible charged 



electron fluid placed in a rigid, neutralizing positive background set to +en . The electron 
density is written as no + n, where n is a small perturbation, and the electronic velocity field 
is denoted v. Finally, we include the electrostatic potential $ and neglect radiation effects. 
The basic equations for the system are the linearized versions of the continuity equation and 
of Euler's and Poisson's equations ^ ]: 
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Here s = (dP / dn) / m = y 3/5 -up is the sound velocity of the liquid. Combining Eqs. 
and introducing the plasma frequency uj p = e 2 n /me we obtain a wave equation for n: 
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For the case of a homogeneous electron liquid confined in a cylinder of radius R we let all fields 
have the dependence f(r,9,z,t) = fi(r) exp[i(W + qz — ut)], with fi being a Bessel function. 
Outside the cylinder $(r) must decay and fulfill Eq. (||) with n = 0, and so $>(r) oc K~i(qr). 
Inside the cylinder, $<(r) can be either decaying, as Ii(nr), or oscillating, as Ji(kr). The lowest 
lying modes are the decaying ones reminiscent of surface plasmons. The oscillation frequency 
uj is found by enforcing the boundary conditions that $(r) and its derivative are continuous at 
r = R, and that the normal component v r of the velocity vanishes at the surface. The solution 
is 

qK[(qR) ,„ _\ „ e 2 n 1D „ , , 2 
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Plasmons of an inhomogeneous electron liquid in a squeezed elliptical cylinder 

Next, to approach the point contact geometry we introduce two perturbations. First, the 
cylinder containing the inhomogeneous electron liquid is squeezed geometrically in a region 
of length 2A around z = 0, i.e. the radius becomes a function of z, say for example R(z) = 
R — SR[1 + cos(7rz/A)]0(A — \z\). Similarly, a static z dependent dip is imposed on the 
positive background charge density n inside the squeezed cylinder, say n (z) = n — 8n[l + 
coa(nz/L)]G(L-\z\). 

In the adiabatic limit where derivatives of R(z) and Uq{z) are neglected, the wave equa- 
tion Eq. (Q) remains separable in cylindrical coordinates, and we make the ansatz n(r, 8, z) = 
Ji(nr)g(z) exp[i(W — uot)], where Ji is a Bessel function and g(z) an arbitrary function to 
be determined. The boundary condition v r (R(z)) = translates into a Neumann bound- 
ary condition J[{nR{z)) = and consequently the "wavenumber" k becomes a function of 
z, k — K n i(z) = j n i/R(z), with j n i being the n'th root of J[{x). Furthermore, uj 2 also 
becomes a function of z, since uJp(z) = e 2 no(z)/me, and similarly for the sound velocity, 
s = s(z) oc nof^) 1 / 3 . As a consequence the wave equation Eq (f|) for n is changed into the 
following eigenfunction equation for g(z): 

- s{z) 2 d 2 z g{z) + [s(z) 2 K 2 nl (z) + uj 2 p (z)} g(z) = uj 2 g(z). (6) 

This is equivalent to Schrodinger's equation (with a position dependent mass) as seen by 
the identifications s 2 <-> Ti 2 /2m and s(z) 2 k 2 i1 (z) + oj 2 (z) <-> V(z). Since ^(0) < uj 2 (±oo) 
bound states, i.e. localized plasmons, may exist. The "effective potential" V(z) is a sum of two 



terms; one, uj^ is dipping down on the length scale L, the other, s 2 fi^, is peaking on the length 
scale A. Depending on the relative strengths, shapes and length scales of the two terms the 
effective potential will appear rather differently. However, for realistic parameters, where the 
density variation dominates, we conclude that localized plasmons may exist in the squeezed, 
inhomogeneous cylindrical electron liquid as shown in Fig. 1. 

The previous considerations dealt with a cylindrical geometry, but it is not difficult to 
approach the 2D case. The trick is simply to use elliptical coordinates (u, v, z) defined by 
(x, y, z) = (77 cosh(w) cos(f ), r\ sinh(w) sin(w), z). The parameter rj relates to the eccentricity of 
the ellipse. With these coordinates the wave equation separates as before. Instead of trigono- 
metric functions of the angle 9 we now obtain the Mathieu functions of the generalized angular 
variable v, and instead of Bessel functions we obtain the modified Mathieu functions of the gen- 
eralized radial coordinate u. By letting the eccentricity 77 tend to infinity we end up with a 2D 
geometry close to the one realized in the quantum point contact experiments. The conclusions 
obtained for the circular cylinder can be restated for the elliptic cylinder, and thus localized 
plasmons are expected to exists in or near the constriction region of quantum point contacts. 



Figure 1: (a) The squeezed elliptical cylinder. High and low densities are represented by dark 
and light shadings respectively, (b) The effective potential V(z) (full line) is determined by the 
parameter values of Ref. [f|. For this potential a solution of the wave equation Eq. (|j) for g(z) 
is found numerically (dashed line). The solution represents a localized plasmon with an energy 
of the order 10 K. 

Plasmon damping 

So far we have treated only the undamped case. In real systems the collective plasmons are 
damped through their interaction with individual electron-hole pairs, the so called Landau 
damping. This effect could be simulated by adding a damping term — v/r to the right hand side 
of the Euler equation Eq. (0) . Instead we will leave the classical level of description and continue 
with a microscopic quantum treatment. The classical level is adequate for demonstrating the 
existence of the collective (almost classical) plasmon excitations, but fails when it comes to 
single particle effects. 

The point contact can be approximated by a ID region, the constriction, connected at each 
end to 2D regions, the contacts. For this 2D-1D-2D model of the point contact we can estimate 
the frequency of the confined plasmon using our insight from the classical calculations: we 
calculate the dispersion relation for an infinite ID-wire and insert the wave vector q c = 2tt/L, 
L being the length of the constriction and hence related to the size of the localized plasmon. 
The long wave length limit of the dispersion relation found by a RPA calculation is 
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Note how the second term under the square root resembles the classical result of Eq. (|^). In 
Rcf. [[J we used this formula successfully to fit the measured activation temperatures mentioned 
in the introduction. 

The confined lD-plasmons will be Landau-damped through their coupling to the 2D-contacts 
outside the constriction. Inserting q c in the RPA expression for the polarizability x 2D we ob- 
tain the following rough estimate of the lifetime t^/uj} of the lD-plasmon coupled to the 
2D-contacts of the 2D-1D-2D model: 

^ « V 2D Im^ D « ^ -^C ^ !" (8) 
The plasmons are seen to be damped, but not over-damped. 
Quasiparticle lifetime 

In the Landauer-Biittiker formalism the conductance is given by single particle properties. Once 
a particle is launched in a given channel of the injecting lead the transmission probability am- 
plitudes are governed by the elastic scattering matrix of the system. For a quasiparticle with a 
finite lifetime it is possible that a particle will decay before completing its traversal of the sys- 
tem. We propose that the observed deviation from perfect quantized conductance is indeed due 
to the finite lifetime of the quasiparticles. Furthermore we suggest that the main contribution 
restricting the lifetime comes from scattering against the localized plasmons. As demonstrated 
above, the localized plasmons provide a well defined finite energy Uojl. Through the Coulomb 
interaction the electrons will scatter against the plasmons and hence the quasiparticle lifetime 
and the transmission properties are affected. The resulting lifetime and additional resistance is 
expected to exhibit an activated behavior, r" 1 oc exp(— T a /T), since a finite energy is needed 
to excite the localized plasmon. We are thus lead to identify the activation temperature with 
the energy of the localized plasmon: T a = hu^/kB- 

Conclusion 

Using a hydrodynamic description of the electron fluid, we have shown that localized plasmons 
with a frequency uj^ are likely to exist near the constriction of a point contact. We have 
sketched how a more complete microscopic quantum calculation may account for a quasiparticle 
lifetime broadening r _1 with a thermal activation behavior t~ 1 oc exp(—hu)L,/kBT). We relate 
this broadening with conductance and are lead to identify the recently measured activation 
temperature T a for conductance deviations with the frequency ujl of the localized plasmon. 
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